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Abstract: 

It has been shown by Gupta and Padmanabhan that the radiation reaction force of the 
Abraham-Lorentz-Dirac equation can be obtained by a coordinate transformation from the 
inertial frame of an accelerating charged particle to that of the laboratory We show that 
the problem may be formulated in a flat space of five dimensions, with five corresponding 
gauge fields in the framework of the classical version of a fully gauge covariant form of 
the Stueckelberg-Feynman-Schwinger covariant mechanics (the zero mode fields of the 
0,1,2,3 components correspond to the Maxwell fields). Without additional constraints, 
the particles and fields are not confined to their mass shells. We show that in the mass-shell 
limit, the generalized Lorentz force obtained by means of the retarded Green's functions for 
the five dimensional field equations provides the classical Abraham-Lorentz-Dirac radiation 
reaction terms (with renormalized mass and charge). We also obtain general coupled 
equations for the orbit and the off-shell dynamical mass during the evolution. The theory 
does not admit radiation if the particle remains identically on-shell. The structure of the 
equations implies that mass-shell deviation is bounded when the external field is removed. 
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Gupta and Padmanabhan 4 have shown that the motion of a charged particle in an 
electromagnetic field can be described in the inertial frame of the particle with a time 
varying non-trivial background metric. Using the general covariant form of the Maxwell 
equations and transforming back to the inertial frame of the laboratory, they obtained 
the Abraham-Lorentz-Dirac radiation reaction term as a consequence of this geometrical 
picture. This result demonstrates that the description of the motion of a charged particle 
in acceleration must include the radiation terms of the Abraham-Lorentz-Dirac equation 5 . 

Alternatively, one can develop the mechanics in a flat space of higher dimension, an 
approach that we shall take. We shall work with the manifestly covariant mechanics of 
Stueckelberg 6 , which provides a description of dynamical systems under the influence of 
forces (which may be represented in terms of potentials or gauge fields) in a framework 
which is Lorentz covariant. This theory admits, on a classical level, deviations from the 
particle's mass shell during interaction, as in quantum field theory. A similar approach 
was used by Mendonga and Oliveira e Silva 7 , who studied the motion of a relativistically 
kicked oscillator in the E, t plane using what they called a "super Hamiltonian." One can, 
in fact derive the relativistic Lorentz force 

= _F M v x v (1) 

from such a Hamiltonian. 

Consider the Hamiltonian 5 ' 6 (we take c = 1 henceforth) 

K= 2M (2) ' 

where x = x^. The Hamilton equations (generalized to the four-dimensional symplectic 
mechanics 6 ) are 

dx>^_dK L _ - eA»(x) 

dp>* _ OK _ dA x {x)p x -eA x {x) 1 ' 

dr dxp, dxp, M 

where r is the absolute (universal) invariant time parametrizing the path of the particle in 
spacetime 6 . Computing -j^- from the first of these, one finds Eq. (1). It moreover follows 
from the first of Eqs.(3) that 

dx^dx^ _ (pV - eA^{x) ) (p M - eA^ (x) ) _ 

dr dr M 2 ' U 

this quantity is absolutely conserved, since K does not depend explicitly on r. It follows, 
since the square of the proper time ds 2 = —dx^dx^, that ds is proportional to dr, inde- 
pendently of the acceleration of the particle. The numerator of (4) is the mass-squared of 
the particle; we infer that this result is associated with the restriction of the particle to a 
sharp mass shell. 
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Taking into account full U(l) gauge invariance, the Stueckelberg-Schrodinger 
equation 6 (including a compensation field for the r-derivative) is 

(il + = (P"-^"fex))fa-eoa, (a: ,x)) ^ M| (5) 

where the gauge fields may depend on r and eo is a dimensionless coupling. The corre- 
sponding classical Hamiltonian then has the form 

K = (P"-eoa^,r))^-e a,(x,r)) _ ^ (g) 

The equations of motion for the field variables are given (for both the classical and quantum 
theories) by 8 

Xd a f a (x,T) = e f(x,r), (7) 

where a, (3 = 0, 1, 2, 3, 5, the last corresponding to the r index, and A, of dimension is 
a factor on the terms f a ^ f a f3 in the Lagrangian associated with (6) (including degrees of 
freedom of the fields) , required by dimensionality, as we shall see below. The field strengths 
are 

fafl = d * a f3 _ d f3 a a^ (8) 

and the current satisfies the conservation law 8 ' 10 

d a j a (x,r) = 0. (9) 
Writing out (9) explicitly (j 5 = p, the density of events in spacetime), 

d 5P + d^ = 0; (10) 

integrating over r on (— oo, oo), and assuming that j 5 (x, r) vanishes 8 at \t\ — > oo, one finds 
that 

d^{x) = 0, 

where (for some dimensionless rf) 

/oo 
drf{x,T). (11) 
-oo 

We identify this J M (x) with the Maxwell conserved current. In ref. 9, for example, this 
expression occurs with 

^(x,r)=^(r)5 4 (x-x(r)), (12) 

and r is identified with the proper time of the particle (an identification which can be 
made for the motion of a free particle). 

Integrating the /(/-components of Eq. (7) over r (assuming / m5 (x, r) — > for r — > ±oo), 
we obtain the Maxwell equations with the Maxwell charge e = eo/rj and the Maxwell fields 
given by 

/oo 
a fI (x,T)dr. (13) 
-oo 
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The Hamiltonian of Stueckelberg 6 and Mondonga and Oliveira e Silva 8 can be recovered 
in the limit of the zero mode of the fields 



/ 



dsa fJ '(x, s)e 



— 1ST 



(14) 



In the zero mode limit, when the Fourier transform of the fields have support only in the 
neighborhood As of s = 0, the vector potential takes on the form a M (x, r) ~ Asd At (x, 0) = 
(As/27rA)A At (x), and we identify e = (As/27rA)eo- The zero mode therefore emerges when 
the inverse correlation length of the field satisfies the relation r/As = 2n\. We remark 
that in this limit, the fifth equation obtained from (7) decouples; the zero mode of the 
r derivative of a^{x,r) vanishes. If the parameter A is independent of the dynamical 
structure of the fields, then the effective width of a M (x, s), when it is well-defined, affects 
the value of the charge e, as well as the relation between the effective Maxwell current and 
the microscopic current j M .This effect, occurring when a Maxwell type theory is a good 
approximation, can be understood as a classical analog of charge renormalization, where 
the effective charge is a function of momentum transfer. 

Again, writing the Hamilton equations for the Hamiltonian (6) , we find the generalized 
Lorentz force 10 



Mi" = eoF v x v + 5 = e (f se i / M v x v + / se // 5 + fext" v x v + f ex S 5) • (15) 



this equation therefore does not necessarily lead to the trivial relation between ds and dr 
discussed above in connection with Eq. (4) . The / M 5 term has the effect of moving the 
particle off-shell. 

In the following we use the Green's functions for (7) to calculate the radiation reaction 
force directly, as, for example, in the derivation of Sokolov and Ternov 11 . In the limit for 
which the particle stays on its mass shell during the interaction, we show that this formula 
reduces to the known Abraham-Lorentz-Dirac formula 9 ' 12 for the Maxwell self-interaction 
problem. We furthermore show that the deviation from mass shell is stable. We shall use 
the retarded Green's function and treat divergences by renormalization of charge and the 
mass parameter M. 

Choosing the generalized Lorentz gauge d a a a = 0, Eq. (7) becomes 



where a = ±1 corresponds to the possible choices of metric for the symmetry 0(4, 1) or 
0(3, 2) of the homogeneous field equations. 

The Green's functions for Eq. (17) can be constructed from the inverse Fourier trans- 
form 



Multiplying this equation by i M , one obtains 



Mx^ = eox^f^ 5 = e (x^f se if^ 5 + x^f ex S 5) ; 



(16) 



Xd a d a a^x,r) = (ad 2 T - d 2 t + V V = -e /(x,r), 



(17) 




(18) 



4 



Integrating this expression over all r gives the Green's function for the standard Maxwell 
field. Assuming that the radiation reaction acts causally in r, we shall restrict our attention 
here to the r-retarded Green's function. In his calculation of the radiation corrections to 
the Lorentz force, Dirac 12 used the difference between advanced and retarded Green's 
functions in order to cancel the singularities that they contain. One can, alternatively 11 , 
use the retarded Green's function and "renormalize" the mass in order to eliminate the 
singularity. In our analysis, we follow the latter procedure. 

The r- retarded Green's function is given by multiplying the principal part of the 
integral (18) by 9(r). Carrying out the integrations (on a complex contour in k; we 
consider the case a = +1 in the following), one finds (this Green's function differs from 
that used in ref. 13, constructed on a complex contour in k°) 



tan 



G (x ,r)= 2 M\ ""(-S^ -^T * 2+r2< ° ; (19) 



(27T) 3 1 1 i I r-y^+^ l -r 2 , ^2 ^ n 

1 2 (r2+x2) | m | r+Vr 2 +x2 | ~ x *( T *+x*) X + r >U. 



With the help of this Green's function, the solutions of Eq. (17) for the self-fields can 
be written, 

a self ^(x,T) = y J d^x'dT'G{x-x',T-T')x^{T')5^{x' -X{r')) 

= y J dr'x^(r')G(x-x(r'),r-r') 
a se i f 5 (x,r) = j J d 4 x'dT'G(x-x',T-r')5 4 (x' - x{r')) 
jdr'Gix-x^r-r') 



(20) 



- ^£ 
~ A 

where we have used (12) (along with j 5 (x, r) = 5 4 (x — x(r))). We have written this Green's 
function as a scalar, acting in the same way on all five components of the source j a ; to 
assure that the resulting field is in Lorentz gauge, however, it should be written as a five 
by five matrix, with the factor — k a kfj/k 2 (/c 5 = k) included in the integrand. Since 
we are computing only the gauge invariant field strengths here, this extra term will not 
influence any of the results. 

From (8) and (15), it then follows that the generalized Lorentz force for the self-action 
(the force of the fields generated by the world line on a point x M (r) of the trajectory) is 

M ^ = ^y dr , (x , '(r)i: I/ (r , )^-i: , '(r)^(r , )9 I/ )G(x-x(r , ))U =a!(T) 

e f j'fau \n<„ . . ( 21 ) 



+ A 



f J dT'(d^-x^r')d T )G(x-x(r'))\ x=x{T) 

+ eo (/ext M vX V + fext^ b) 

We define u = (x„(r) — x M (r'))(x /i (r) — x M (r')), so that 



a M = 2(x M (r)-x M (r'))^. (22) 
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Eq. (21) then becomes 

MF = 2y y dr'{x u {T)x v {T'){x^{T) - x"(t')) 

- x»(t)x»(t')(x„(t) - ^(r'))}^^* - x(/),r - t')\ x=x{t) ) 

+ | | dr / {2( a ^(r)- a ^(/))| ; -^(r / )a r }G(x-x(/),r-/)| jB=a!(T) ). 

+ eo (/ex/ ^ + fext 1 5) 

We now expand the integrands in Taylor series around the most singular point r = r' . 
In this neighborhood, keeping the lowest order terms in r" = t — t', the variable u reduces 
to u = x^x^r" 2 . We shall also use the following definition; 

e=l + x"x^ (24) 

a quantity that vanishes on the mass shell of the particle (as we have pointed out above). 
In this case the derivatives of (19) take the form 



dG „ 0(T")A(e) 



du (27r) 3 r //5 
dG ^ 9(T")f 2 (e) <*(r")/ 3 (e) 
dr" (2tt)V' 4 (2tt)V 3 



(25) 



where we have used the following definitions: 

e < 0: 



3tan-V5£) _ 3 2 
(-e)i e 2 (l-£) e(l-e) 2 

, , N Stan-^v^i) 1 2-e / N 

/2(£) = — ? " ^(T^T) ( 26fl ) 

/3(e) = 3 + 



(_ £ )2 e(l-e) 



£ > 0: 



3 ln |l±V£| 

/i(e)= 2 ^' -^7^ + 



tln|i±^l 



/a (e) = l^^l - 1 - ^-L, (266) 
\Hr^r\ 1 

/3(e) = 3^ + 
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For either sign of e, when e ~ 0, 

~ g + ye, 

/ 2 (e)~-|-^e, (26c) 

2 ^ 4 

/3(e) ~ 3 + 

One sees that the derivatives in (25) have no singularity in e at e = 0. 
From (8) and (20), we have 

fself 1 " 5 (x(t),t) = 

e J dr'{2(x^r)-x^r'))^--x^r')d T }G(x-x(r'),r-r')\ x=x{T) , 



(27) 



We see (from (25)) that the main contributions to the integrals come from small r" . We 
may therefore expand x M (r) — x M (r') and x^(t) — x M (r') in (27) in power series in t", and 
write the integrals formally with infinite limits. 

Substituting (27) into (16), we obtain (note that x M and its derivatives are evaluated 
at the point r, and are not subject to the r" integration), after integrating by parts using 
S(t") = ^0(t"), 

(/1-/2-3/3). .„ (J/l - /2 " 2/3) . .. v 



. r _ J4_ r- (/l -/ 2 -3/3) (j/i-/ 

' A(2;;)-\/ /' 1 t"4 ^ (2g) 



^ a^z }#(r ) +e x^f ext tl 5 . 



1 T //2 

The integrals are divergent at the lower bound r" = imposed by the ^-function; we 
therefore take these integrals to a cut-off ji > 0. Eq.(28) then becomes 

M = H f (/i-/ 2 -3/ 3 ) (^1-/2-2/3) 
2 A(2tt)^ 3^ ^ * V (29) 

XjyX I + eoX^fext^ 5. 



Following a similar procedure, we obtain from (23) 



(/1 - /2 - 3/3) - / 2 - 2/3) (|/i - \ f 2 - \h) ■ 



+ 



+ e {f ex S„X» + f ex / 5 ). 



2^ 2 



2^ + 



(30) 
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Using (29) to substitute for the coefficient of the -\ term in (30) , we obtain (for e 7^ 1) 



2 



ufext 5 . /" u • v 1 /■ u 

H : h e Jext vX + e Jext 5, 

1 — £ 



(31) 



where 

F(e) = £(l-e) + (I/i - I/ 2 - I/ 3 ). (32) 

Here, the coefficients of x^ have been grouped into a renormalized (off-shell) mass 
term, defined (as in the procedure of Sokolov and Ternov 11 ) as 

M(e) = M + ^- + \h ~h- 2/3] (33) 



where, as we shall see below, 



A(2tt)V [6V 



can be identified with the Maxwell charge by studying the on-shell limit. 
We now obtain, from (31), 

M(e)x» = -l^^-ex^ + F(e)e 2 Vx fl + -^— Xv x v x»\ 
2 1 — e 1 — e 



(35) 



We remark that when one multiplies this equation by x^, it becomes an identity (all of 

the terms except for eof ext ^x u may be grouped to be proportional to {^zf + S^j); one 

must use Eq. (29) to compute the off-shell mass shift e corresponding to the longitudinal 
degree of freedom in the direction of the four velocity of the particle. Eq. (35) determines 
the motion orthogonal to the four velocity. Equations (29) and (35) are the fundamental 
dynamical equations governing the off-shell orbit. 

We now show that the standard relativistic Lorentz force, with radiation corrections, 
can be obtained from these equations when \xe « e « 1 and e and fj? xt 5 are small. In 
this case, Eq. (29) becomes 

(M — ) - = e 2 {-^ + —x v x") (36) 

V 15//' 2 I 15^ 2 15 J y ' 

The left hand side can be neglected if 



.e 2 



M/(-) (lie) « e. (37) 
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We shall see below that we must have 0.68e 2 //j < M for stability of e, but if e 2 / is 
not too small, the inequality (37) is consistent with our assumed inequalities, and it then 
follows that 

4e//j 2 = x v x '. (38) 

If, furthermore, e is small, then 

x jj, x — s, — i (^39^) 

the known expression associated with radiation. Since e/fx 2 may be appreciable even if 
e is small, the inequalitites we have assumed can admit a significant contribution of this 
type. Under these conditions equation (34) becomes, 

M rm x» = ^e 2 {i' M - x v x v x»} + e f ex S „x v , (40) 

where M ren = M(e)| e=0 = M + e 2 /3/i. 

This result is of the form of the standard relativistic Lorentz force with radiation 
reaction. 9 ' 11 ' 12 ' 14 

We now study the stability of the variations of the off-shell parameter e when the 
external field is removed. First, we construct an equation of motion for s. We define the 
functions 

F 1 (e) = ^(e-l)(f 1 -f 2 -3fs) 

F*{e) = i-(i/ 1 -/ 2 -2/ 3 ) (41) 

Fz( £ ) = g/i _ 2^ 2 ~~ 2^ 3 
so equation (29), in the absence of external fields, becomes: 

= e 2 {F 1 (e) + F 2 (e)e + F 3 (e)x^}. (42) 

Solving for the explicit x derivatives in (42) and differentiating with respect to r, one 
obtains 

x"x> + =y 3 {F^e 2 + + F 2 ) - F[e} 



Together with 

one finds, from (43), 



Fo ( M , 
^{^ + 2 ? e-F 1 }e = J ff. 



. A 4 , O" "■A 1 _ 1 ••• 

X^X -\- oX^X — — £ 



(43) 



V M = \t ~ \H(£,£,£) (44) 



Multiplying Eq.(35) by (with no external fields) and using (42) and (44), we obtain 

e - A(e)e + B(e)e 2 + C(e)e - D(e) = 0, (45) 



where 



2 , M . 2M(e) 
^ e ) = -£r(o^ + ^) + 



F 3 v 2e 2 A} e 2 F( £ )' 
2F'M , 2FL, 2 1 ,M , M(e) 



Ff v 2e 2 ' F 3 l- £ F 3 v 2e 2 J e 2 F(e)l-e' 
4M(g) 1 ( M 2 

" (l-e)F 3 + T 3 15 



D = 



4M(e) Fi 
e 2 F(£)Fs"' 



We first study the possibilty of having a solution of the form e = Eq, a constant. In 
this case e = implies, 

Since all the derivatives of £ are zero we also find from (44) , 

a^cc =0 

. Multiplying eq.(35) by x M and substituting these last two results we get 

liy l^iiy lis tAJ ^ \J 



From (42) we find then that 



x„a; = — — = 0. 
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From (41) and (266), one sees that this equation can be satisfied only if e = (Fi = 0) or 
e = l(F 3 = oo). 

Since e = we find that i 2 i 2 = |x| 2 |x| 2 cos 2 #. Together with x^ = i 2 - |x 2 | = this 
implies 

i 2 p = \ x \ 2 i 2 cos 2 e 

The solution i = \x\ =0 implies = const . The other solution i 2 = |x| 2 cos 2 6> implies 
that |x| 2 (l — cos 2 6) = e — 1; since the left hand side is positive, e cannot be zero, and the 
only possibility for a constant solution is then e = 1, motion on the light cone. We shall 
show below that the trajectory cannot reach this bounday. 

The mass shell condition e = 0, in the theoretical framework we have given here, 
implies that the particle motion must be with constant velocity, and that no radiation 
(x fX x fl = 0) is possible, i.e., in order to radiate, the particle must be off-shell. This result 
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is also true in the presence of an external field. In particular, it follows from Eq. (29) that 
for e = 0, 

2 ... 2 

15// ^"15^ (47) 

From Eq.(15), however, it follows (in case s = 0) that x M / e xt M 5 = —i-ufseif^b, so that 
the nonvanishing value of i^f corresponds only to a self-acting field f se if^ 5 (driven by 
fext 1 5), and not to radiation. 

We now show that, in general, e is bounded when the external fields are turned off. 
For the case e < the function F 3 is zero at e = —0.735. In this case eq.(42) becomes 

If e > at this value of e, then £ cannot cross this boundary. Since Fi(— 0.735) = , 
F 2 (-0.375) = -^p, this condition implies that 

e 2 

A« > 0.68—. 

r M 

Setting M, e equal to the electron mass (the lowest mass charged particle) and charge one 
finds that \i > 10 cut-off of reasonable size for a classical theory. 

We now show that e is bounded from above by unity. The full classical Hamiltonian, 
obtained by adding the contribution of the fields to the expression on the right hand side 
of (6), is a conserved quantity. In the absence of external fields, all the field quantities are 
related to the source particle through the Green's functions. In the absence of external 
fields, as the particle motion approaches the light cone, there are infinite contributions 
arising from the fields evaluated on the particle trajectory. In this case, it follows from 

(4) that (P^-^(x,r))^-e a^x,r)) = Q ^ ^ gelf _ field term ig ^ singular than the 

fnvf^ and ^5/^5 terms, which involve derivatives of the Green's functions, as in (25), 
squared. As seen from (266), the most singular contribution arises from fi 2 . Since the 
total Hamiltonian K is conserved, the coefficient of this singularity must vanish. The 
coefficients involve just e (and its square) and e; one finds a simple nonlinear differential 
equation for which only e = can be a solution. It follows that the conservation law 
restricts the evolution of e to values less than unity, i.e. the particle trajectory cannot 
pass through the light cone. 

This bound manifests itself in the structure of the differential equation (45) for e. In 
the limit that e — > 1, the coefficients A, B, C, D are all finite; however the behavior of the 
linearized solution depends on the derivates of these coefficients, and, in this limit, B' is 
singular, driving the solution away from the light cone. 

Numerical studies are under way to follow the motion of this highly nonlinear system 
both in the presence and absence of external fields. 
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